Correlations between large prime numbers 
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It is shown that short-range correlations between large prime numbers (~ 10^ and larger) 
have a Poissonian nature. Correlation length C, — 4.5 for the primes ~ 10"" and it is increasing 
logarithmically according to the prime number theorem. For moderate prime numbers (~ 10'') the 
Poissonian distribution is not applicable (while the correlation length C, surprisingly continues to 
follow to the logarithmical law. A chaotic (deterministic) hypothesis has been suggested to explain 
the moderate prime numbers apparent randomness. 
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Introduction. The prime number distribution is ap- 
parently random. The apparent randomness, hovifever, 
can be stochastic or chaotic (deterministic). The na- 
ture of the randomness is still unsolved problem. Using 
probabilistic methods we should take into account that 
there are different levels of averaging. One can consider 
the probabilities as representing a first level of averaging 
while the correlation functions represent a second one. 
Although the first level contains more details the second 
one can provide a more robust picture. In order to use the 
second level of averaging let us define a binary function 
v{n) of integers n = 3,4,5,...., which takes two values 
or -1 and changes its sign passing any prime num- 
ber. This function contains full information about the 
prime numbers distribution. Due to the prime number 
theorem, which states that the "average length" (r^) of 
the gap between a prime p and the next prime number is 
proportional (asymptotically) to Inp, the v{n) is a statis- 
tically non-stationary function. However, for sufficiently 
large n this function can be considered (in a certain level 
of averaging) as a stationary function in the windows cen- 
tered at no with width A/riQ <C 1. An heuristic estimate 
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FIG. 1: Autocorrelation of the v{n) function in a window with 
the width A = lO** centered at no ~ 10^. 
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FIG. 2: Autocorrelation of the v{n) function in a window with 
the width A = 10^ centered at no ~ 10^. 
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for A/p ^ 1 provides a basis for such consideration. 
We will show (numerically) that on the level of the 
correlation function this idea indeed seems to be rather 
well applicable for sufficiently large prime numbers. On 
the level of the probabilities, however, application of 
this idea demands certain additional average (after that 
the results obtained by the two different methods can be 
compared). 

Large prime numbers. Before performing this analy- 
sis let us recall certain simple properties of the random 
telegraph signal. The random telegraph signal is a bi- 
nary process x{n) which takes two values -f 1 or -1, and 
has Poissonian distribution of the values of n where x{n) 
changes its sign. For the stationary case the Poisson dis- 
tribution results in the exponential distribution of gaps 
T between consecutive values of n where x{n) changes its 
sign 
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FIG. 3: Distribution of the gaps between consecutive primes 
for v{n) in a 'window of stationarity' with the width A = 10* 
centered at no ^ 10^. 
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FIG. 4: As in Fig. 3 but in a window with the width A = 10^ 
centered at no ~ lO'^. 

with constant value of > 0. The autoeorrelation func- 
tion of the statistically stationary random telegraph sig- 
nal x{n) 

C{1) = {x{n)x{n + I)) - {x{n}){x{n + I)) = e'^'^ (3) 

where the angular brackets denote averaging over real- 
izations and 



Figure 1 shows a correlation function computed for 
the v{n) function in a 'window of stationarity' with the 
width A = 10^ centered at no ~ 10^. The semi-log 
axes has been used in this figure in order to indicate 
exponential decay, Eq. (3), of the autocorrelation 
function (the straight line). The correlation length 
C ~ 4.5. Figure 2 shows analogous autocorrelation 
function but in a window with the width A = 10^ 
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FIG. 5: Autocorrelation of the v{n) function in a window with 
the width A = 4000 centered at no ~ 18000. 

centered at no ^ 10^. The correlation length C, ~ 6.3 
in this case. Figures 3 and 4 show the distribution of 
the gaps T between consecutive primes: P(r), for v{n) 
in a 'window of stationarity' with the width A — 10^ 
centered at no ~ 10^ and in a window with the width 
A = 10® centered at no ~ 10^ respectively. The straight 
lines in these figures represent the best fit to the data 
(also certain kind of averaging) and the slopes of these 
best fits provide us with the exponents of the Eq. (2): 
6 ~ 9 and ~ 12.8 respectively. Now one can compare 
the values of the correlation length C with values of the 
and with the Poissonian relation Eq. (4). One can 
see that in the second level of averaging correspondence 
to the Poissonian distribution of the primes in the 
considered windows is rather good. Similarly to the 
" average gap" (see above) the second level's parameters 
C and © = 2C are proportional to In no due to the prime 
number theorem. 

Moderate prime numbers. Figure 5 shows an expo- 
nentially decaying part of the autocorrelation computed 
for the v{n) function in a 'window of stationarity' with 
the width A = 4000 centered at no ~ 18000 (moderate 
values of the primes). Naturally, the exponential range 
here is rather short. However, the computed value of the 
correlation length Q ~ 3.83 follows quite precisely to the 
logarithmic dependence on no. It is surprising if we take 
into account that the distribution of the gaps between 
consecutive primes (the parameter 0) does not, as one 
can see from figure 6 (the ~ 9 here). This means, 
in particular, that the v{n) function in this window is 
not a random telegraph signal. Something wrong is go- 
ing with the Poissonian hypothesis for these values of the 
primes. But if previously we related the logarithmical de- 
pendence of the correlation length Q on no just with the 
logarithmical dependence of 0, then how can we explain 
it now? 

The problem with the Poissonian distribution for the 
moderate primes may be related to the strong and persis- 
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FIG. 6: Distribution of the gaps between consecutive primes 
for v{n) in a 'window of stationarity' with the width A — 4000 
centered at no ~ 18000. 
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FIG. 7: Spectrum of the oscillations shown in Fig. 6 (after a 
detrending) . 



tent periodic oscillations in the distribution of the gaps. 
Figure 7 shows spectrum of these oscillations (seen in 
Fig. 6) after a detrending, and the strong peak in Fig. 
7 corresponds to the period equal to 6 (as well as for 
the larger no). In the Refs. [H-@ this period has been 
studied both numerically and theoretically. We are more 
interested, however, in nature of the randomness of the 
moderate primes. The spectrum of the function v{n) in 
the moderate window of stationarity, shown in figure 8 
in the semi-log scales, may help in solving this problem. 
One can see a pronounced range of an exponential decay 
of the spectrum. Many of the well known chaotic attrac- 
tors ('Lorenz', 'Rossler', etc.) exhibit the exponentially 
decaying spectra and this spectrum is considered as a 
strong indication of a chaotic behavior [3|-0- Here, for 
comparison with the Fig. 8, we will consider a chaotic 
spectrum generated by the Kaplan- Yorke map Q . In the 
Langevin approach to Brownian motion the equation of 
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FIG. 8: Spectrum of the function v{n) in a 'window of sta- 
tionarity' with the width A — 4000 centered at no ~ 18000. 

motion is 

y = —fy + F{t) (5) 

where the fluctuating kick force on the particle is a Gaus- 
sian white noise: F{t) = J^nVnSU — nr) and y{t) to take 
values in R"^. One can assume Q that the evolution of 
the kick strengths is determined by a discrete time dy- 
namical system T on the phase space and projected onto 
i?™ by a function g: 

Vn = 5(a:„_i), Xn+l = TXn (6) 

Then the solution of Eq. (5) is 

y[i) = e-^(*-"^)y„ (7) 

where n equals the integer part of the relation t/r and 
the recurrence 

Xn+i = Tx„, yn+1 = ayn + gixn) (8) 

provides value of yn (with a = e~^'^). In certain sense 
the dynamical system (8) is equivalent to the stochastic 
differential equation (5). In the generalization related to 
the Eq. (8) the force F{t) can be considered as a non- 
Gaussian process which is determined by g and T. 

The Kaplan- Yorke map is a particular sim- 

ple case for this generalization: 

Tx ~ 2x {mod 1), g{x) = cos47ra; (9) 

An optimal computational algorithm for the Kaplan- 
Yorke map is 

a„+i = 2a„ modip), Xn+i = ttn/p, yn+i = a2/„+cos47r2: 

(10) 

where p is a large prime number. 

Figure 9 shows spectrum of a chaotic solution of the 
Kaplan- Yorke map {a = 0.2). We used the semi- log axes 
in this figure in order to indicate exponential decay of 
the spectrum (the straight line). 
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FIG. 9: Spectrum of a chaotic solution of the Kaplan- Yorke 
map. 



One can compare Fig. 9 and Fig. 8 in order to conclude 



on the chaotic nature of the "window defined" 
tion of the moderate primes. 
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A rigorously minded reader can find a lot of good math- 
ematics related the above discussed approach: such as 
Gallagher theorem for short intervals 12] (about Poisson 
distribution in the limit case of logarithmically short win- 
dows), Cramer's probabilistic model [13] and the difficul- 
ties of this model in the short intervals (Maier's theorem 
H'ISl'IS)- Since present consideration was devoted 
to the correlations between large primes and the corre- 
lation length ( is rather small (see Figs. 1,2 and 5) the 
consideration was restricted by the windows of statistic 
stationarity where the Poisson distribution seems to be 
reigning on the second level of statistical averaging at suf- 
ficiently large primes. Above these windows distributions 



of the primes can be different (see, for instance, Ref. [17] 
and references therein). For moderate values of primes 
the Poisson distribution seems to be not applicable and 
a chaotic hypothesis has been suggested instead. 



]1] M. Wolf, Applications of statistical mechanics in number 

theory, Physica A, 274, 149-157 (1999) [9] 

[2] P. Kumar, P.Ch. Ivanov, and H.E. Stanley, Infor- 
mation Entropy and Correlations in Prime Num- JlO] 
hers, |arXiv:cond-mat / 03031 10] ;a4 [cond-mat.stat-mech], 
(2003). 

]3] S. Aresa, and M. Castro, Hidden structure in the ran- Jll] 
domness of the prime number sequence?, Physica A, 360, 
285-296 (2006). J12] 

]4] D.E. Sigeti, Survival of deterministic dynamics in the 

presence of noise and the exponential decay of power J13] 
spectra at high frequency, Phys. Rev. E, 52, 2443; Phys- 
ica D, 82, 136 (1995). 

]5] N. Ohtomo, K. Tokiwano, Y. Tanaka et. al.. Exponential J14] 
Characteristics of Power Spectral Densities Caused by 
Chaotic Phenomena, J. Phys. Soc. Jpn., 64 1104 (1995). J15] 

]6] J. D. Farmer, Chaotic Attractors of an Infinite- 
Dimensional Dynamical System, Physica D, 4, 366 J16] 
(1982). 

J7] U. Frisch and R. Morf, Intermittency in nonlinear dy- 
namics and singularities at complex times, Phys. Rev., J17] 
23, 2673 (1981). 

[8] J.L. Kaplan and J. A. Yorke, in Functional Differential 
Equations and Approximations of Fixed Points, Lecture 
Notes in Mathematics, 730, p. 204. (New York: Springer 



1979). 

C. Beck, Ergodic Properties of a Kicked Damped Parti- 
cle, Commun. Math. Phys., 130, 51 (1990). 
R.V. Jensen and C.R. Oberman, Calculation of the Sta- 
tistical Properties of Strange Attractors, Phys. Rev. 
Lett., 46, 1547 (1981). 

J. C. Sprott, Chaos and Time-Series Analysis (Oxford 
Univ. Press, 2003) 

P.X. Gallagher, On the distribution of primes in short 
intervals, Mathematika, 23, 4-9, (1976). 
H.Cramer, On the order of magnitude of the difference 
between consecutive prime numbers. Acta Arith., 2, 23- 
46 (1936). 

H. Maier, Primes in short intervals, The Michigan Math. 
Journal, 32 221 (1985). 

A. Granville, Proc. Intern. Congress of Math., I, 388 
(1995) (Zurich, Switzerland, 1994). 

H.L. Montgomery, and K. Soundararajan, Primes in 
short intervals, Commun. Math. Phys., 252, 589-617 
(2004). 

K. Soundararajan, The distribution of the prime num- 
bers, NATO Science Series II: Mathematics, Physics and 
Chemistry, , 237, 59-83 (2007). 



